arXiv: 1505.03409v2 [physics.data-an] 20 Jan 2016 


Search for a Lorentz invariant velocity 
distribution of a relativistic gas 

Evaldo M. F. Curado^’^, Felipe T. L. German!^ 
and Ivano Damiao Soares^ 

^ Centro Brasileiro de Pesquisas Fisicas - CBPF 
^ National Institute of Science and Technology 
for Complex Systems 
Rio de Janeiro, Brazil 


Abstract 

We examine the problem of the relativistic velocity distribution 
in a 1-dim relativistic gas in thermal equilibrium. We use numeri¬ 
cal simulations of the relativistic molecular dynamics for a gas with 
two components, light and heavy particles. However in order to ob¬ 
tain the numerical data our treatment distinguishes two approaches 
in the construction of the histograms for the same relativistic molec¬ 
ular dynamic simulations. The hrst, largely considered in the liter¬ 
ature, consists in constructing histograms with constant bins in the 
velocity variable and the second consists in constructing histograms 
with constant bins in the rapidity variable which yields Lorentz in¬ 
variant histograms, contrary to the first approach. For histograms 
with constant bins in the velocity variable the numerical data are fit¬ 
ted accurately by the Jiittner distribution which is also not Lorentz 
invariant. On the other hand, the numerical data obtained from his¬ 
tograms constructed with constant bins in the rapidity variable, which 
are Lorentz invariant, are accurately fitted by a Lorentz invariant dis¬ 
tribution whose derivation is discussed in this paper. The histograms 
thus constructed are not fitted by the Jiitter distribution (as they 
should not). Our derivation is based on the special theory of relativ¬ 
ity, the central limit theorem and the Lobachevsky structure of the 
velocity space of the theory, where the rapidity variable plays a crucial 
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role. For jc? <C 1 and 1//3 = kBT/rriQC^ <C 1 the distribution tends 
to the Maxwell-Boltzmann distribution. 


1 Introduction and summary of results 

In physics, it is difficult to overestimate adequately the importance 
of the Maxwell-Boltzmann (MB) distribution of velocities for gases, 
introduced by Maxwell in 1860 [T]. It was the first time that a proba¬ 
bility concept was introduced in a physical theory, as the existing the¬ 
ories at the time were purely deterministic like Newtonian mechanics 
and wave theory. Actually the work of Maxwell was the starting point 
for Boltzmann to elaborate his research program on the evolution of a 
time dependent distribution of velocities for gases, culminating in the 
articles of 1872 [2] and 1877[3], among other important papers, which 
are amid the fundamental cornerstones of the modern kinetic theory 
of gases and of statistical mechanics. 

With the implicit use of the atomic theory of matter (at that time 
a controversial theory), the new concept of entropy was established 
having also as its starting point the introduction, by Maxwell, of the 
concept of probability. Since then the Maxwell-Boltzmann (MB) dis¬ 
tribution played a fundamental role in the statistical description of 
gaseous systems with a large number of constituents. Actually, in 
many cases it is considerably simpler, and even as accurate as, to 
use the MB distribution instead of Bose-Einstein and Fermi-Dirac 
distributions [1]. However a clear limitation of the MB distribution is 
its nonrelativistic character, encompassing velocities larger than the 
velocity of light in contradiction with the special theory of relativity. 

In the present paper we introduce a 1- dim Lorentz invariant (LI) 
distribution of velocities for a relativistic gas in thermal equilibrium 
which has the MB distribution as a limit for velocities much smaller 
than the velocity of light (with correspondingly relatively small tem¬ 
peratures). A discussion about this distribution in the 1-dim and 
3-dim cases was given by two of us in [5]. Our derivation was based 
on three pillars: the special theory of relativity, the central limit the¬ 
orem and the fact that the velocity space of special relativity is a 
Lobachevsky space, where the additivity of velocities in the Galilean 
relativity is transferred to the additivity of rapidities. 

We have made a numerical simulation of the relativistic dynamics 
for a 1-dim relativistic molecular gas with two components (light and 
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Figure 1: Plot of the equilibrium velocity distributions of heavy particles 
{N 2 = 2, 500) obtained from the same numerical simulation for the relativis¬ 
tic molecular dynamics of a 1-dim gas with N = 5,000 particles: {left) the 
points are obtained from a histogram of constant velocity bins, which is not 
LI, and the continuous curve is the best fit of the points to the Jiittner dis¬ 
tribution do]), with best fit parameter mg/Sj — 0.8570 (rms error ~ 0.017). 
{right) The points are obtained from a histogram of constant rapidity bins, 
which is LI, transformed into a normalized histogram in the velocities (via 
eqs. fITU]) and fITT]) ) with non constant bins but still LI. The continuous curve 
is the best £t of the points to the LI velocity distribution (|8]) with best fit 
parameter f3 ~ 1.4070 (rms error ~ 0.004). 
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heavy particles) in the same vein of Cubero et al. [6]. We used Ni = 
2, 500 light particles of mass rriQ and N 2 = 2,500 heavy particles of 
mass 2mo, with a total of 100 relativistic simulations (taking the box 
as the frame of reference in the simulations), and collected the data 
fixing the number 60 of points to be obtained in the histograms. The 
initial random distributions of the velocity of the particles must satisfy 
— 1 < Vi/c < 1 with Vi = 0 (so that the box has zero velocity). 
The control of how much relativistic is the system is made by how the 
interval of the initial random distributions is closer to 1: the random 
choice of \vi/c\ in the interval [0.75,1) is considered less relativistic 
than in the interval [0.85,1). The four-momentum Pa of the particles is 
conserved for the whole relativistic evolution, with the relative speed of 
the colliding particles being a constant in the collision. The numerical 
methods follow closely [TIEIE]. 

To proceed, for the same relativistic simulations we adopted two 
distinct approaches in the construction of the histograms of the distri¬ 
bution of velocities in the thermalization limit: the first, histograms 
with constant velocity bins (which are not LI) and the second, his¬ 
tograms with constant rapidity bins (which are LI). In the second ap¬ 
proach, from the constant rapidity histograms we obtain histograms 
of the velocities by using the inverse transformations of Eqs. m and 
m- These transformations change the relative size of the bins, with 
the bins no longer equally spaced in a velocity scale as they were in a 
rapidity scale. However these histograms still remain LI. 

It is crucial to remark here that the use of the rapidity is not merely 
a change of variables. The core of our new approach in this paper is 
that we treat the numerical data obtained from the molecular dy¬ 
namics simulations by constructing the histograms with constant bins 
in the rapidity which is a LI procedure. This yields a LI probabil¬ 
ity distribution that has a completely distinct measure from previous 
treatments in the literature (cf. [6]) where the construction of his¬ 
tograms is made with constant bins in the velocity, a procedure that 
is clearly not LI and results in a completely distinct distribution. The 
distribution made with constant bins in the velocity cannot be trans¬ 
formed into the LI distribution obtained in this paper by a change of 
variables. 

Part of our main results are displayed in Figs. [1] and [2] where we 
show the equilibrium velocity distributions derived from our numerical 
simulations for the light and heavy components of the 1-dim relativis¬ 
tic gas, respectively. For the same simulation, the points are obtained 
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Figure 2: Plot of the equilibrium velocity distributions of light particles 
{Ni = 2, 500) corresponding to the same numerical simulation of Figs. [U 
{left) the points are obtained from a histogram with constant velocity bins, 
which is not LI, and the continuous curve is the best £t of the points to 
the Jiittner distribution ([H]), with best fit parameter m^fij ~ 0.8862 (rms 
error ~ 0.009). {right) The points are obtained from a histogram of constant 
rapidity bins, which is LI, transformed into a normalized histogram in the 
velocities (via eqs. ca and ffTT]) i with non constant bins but still LI. The 
continuous curve is the best fit of the points to the Lorentz invariant velocity 
distribution ([H]) with best fit parameter ft ~ 0.9254 (rms error ~ 0.008). 


either from a histogram of constant velocity bins (left fignres), which 
is not LI, or from a histogram of constant rapidity bins (right fignres) 
which is LI. 

In Figs. [T]and[2] (/e/t), where histograms of constant velocity bins 
are adopted, the continuous curves correspond to the best fit of the 
Jiittner distribution (l9|) [TO]. The accurate fits reproduce a result of 
Cubero et al. [6| which is a standard reference in the literature on 
the subject, being considered as a clear evidence that Jiittner’s is the 
correct relativistic distribution although not being a LI distribution. 
With the same procedure this result was later also verified for the case 
of 2-dim [TTl [12] and 3-dim [13] distributions. 

In Figs. [T]and| 2 |(rz 5 /it), where histograms of constant rapidity bins 
are adopted (properly transformed to a histogram where the abscissa 
is the difference of velocities), the continuous curves correspond to the 
best fit of the LI velocity distribution (|8|l. 

These two contrasting results, derived from the same numerical 
simulations, are connected to the distinct particular choices made on 
constructing histograms, namely, whether we choose the histogram to 
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Figure 3: Plot of the velocity distribution of Figs. [2] (right): the dashed 
curve is the best £t of the points to the Jiittner distribution ([9]), with best 
£t parameter mg/dj — 2.0406 (rms error ~ 0.04). We see that while the LI 
velocity distribution fits quite accurately the points, the Jiittner distribution 
cannot be a candidate for a LI distribution. 


be LI or not, a fact that has not yet been taken into account in the 
literature. 

On the other hand the graphs of Fig. [3l done for light particles, 
show that the points obtained from a histogram of constant rapidity 
bins (which is LI) fit quite well the LI distribution ([8|) (continuous 
line) and give evidence that the Jiittner distribution ([9|) cannot be a 
candidate for a LI velocity distribution (dashed line). 


2 A proposal for a 1-dim LI velocity 
distribution 

We start by presenting a derivation of the MB velocity distribution 
which differs from the derivation used by Maxwell but which will be 
more appropriate in our derivation of a LI velocity distribution. Let us 
first consider the addition of velocities in the Galilean space. We know 
that the velocities add according the rule v = Vj . Assuming that 
the velocities of the particles Vj are random variables, with zero mean, 
and considering that the sum is over a very large number of particles, 
we have - by the central limit theorem m - that the probability dis¬ 
tribution of velocities for the random variable u = (1/y/N) v, ap¬ 
proaches the distribution P(i') oc exp(—6z^^) or P(v) oc exp(— ftv^/V), 
recovering thus the famous MB distribution if we plot vTV P(v) versus 
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v/VN. 

As well known, in the special theory of relativity the relative veloc¬ 
ity V of two particles with arbitrary velocities vi and V 2 , with respect 
to a hxed inertial frame, is given by 


V = 


Vi - V2 + (7('^^2) - l)(v2/'fi)[vi • V2 - vl] 
7(?;2)(1 - Vi • V2 /c2) 


( 1 ) 


The above expression also holds with the interchange vi V 2 . The 
square of the modulus of the relative velocity is given by 


^2 _ - V2)^ - (l/c2)[vi A V 2 ]" 


(1 - Vi • V 2 /C 


. 2'\2 


( 2 ) 


which is symmetric with respect to vi and V 2 . It is also a well 
known fact that the square of the relative velocity ([2|) is invariant 
under Lorentz transformations [15]. Let us now consider a fixed inertial 
reference frame, say the laboratory frame. For simplicity, we will 
assume one dimensional only. Let the velocity of two particles be vi 
and U 2 , as measured in this inertial frame. The relative velocity of the 
two particles is given from ([2]) as 


V = 


Vl + V2 

1 -I- vxv^jc? 


(3) 


which is invariant under Lorentz transformations. However m can 
be rewritten as 


l-|-u/c / 1 -|-Ul/c\ / 1 -|-U 2 /C 


1 — v/c V 1 — Ul/c/ V 1 — U 2 /C/’ 
and can be extended to any number of particles, 
1-l-u/c Yj/l + Vi/c 


n(f 


l — v/c \ 1 — Uj/c/ 

' I ' 

Taking the logarithm on both sides of Eq. (|5|) and defining 


(4) 


(5) 


1 (1 + Vi/c 

(Tj = - In -— 

2 \l-v^/c 

we can express (l5|) as 


= tanh ^{vi/c) , Oi G (— 00 , 00 ), 


a = 




( 6 ) 


(7) 


7 











The variable defined in Q is the rapidity associated with the velocity 
Vi and will play a fundamental role in the remaining of the paper. 
Accordingly the rapidity of relatives velocities are therefore quantities 
which are LI and are additive in the arithmetic sense. 

Let us then consider a relativistic gas with a large number of par¬ 
ticles, each with a velocity ^ as measured with respect to the iner¬ 
tial frame of the laboratorjo. We concur that any Boltzmann-type 
equation (relativistic or not) that give rise to a universal stationary 
velocity distribution implicitly assumes the presence of a spatial con¬ 
finement, thus singling out a preferred frame (cf. Cubero et al.[ 6 ] 
and references therein). However as we will see the LI distribution 
function derived here will be independent of such a singling out of 
frame, as the variables involved are LI. If we assume that the veloc¬ 
ities Vi are independent and random variables, with zero mean, the 
variables Ui’s are also independent and random with zero mean, and 
we have - in accordance with the central limit theorem - that the 
probability distribution for the variable s = a/y/N in an interval s 
and s + ds approaches P{s)ds = Cie~^^^ds or, denoting (5 = /3/A^, 
P{a)da = da/y/N, where Ci = y^/S/vr is a normalization 

constant. Using ([ 6 ]) and that da = 'y'^dvjc, where 7 = l/y^l — 
we can write the probability distribution for the velocities of a 1 -dim 
relativistic gas as 


P{v)dv = Cie 




tanh ^( h / c ) 


j'^dv/cy/N. 


( 8 ) 


The factor in ([ 8 ]) is the LI 1-dim line element. The adimensional 
parameter /3 = ijiqc?/kbT, where kb is the Boltzmann constant, de¬ 
creases with the temperature T. Such temperature is obviously LI. 

The distribution Q is a 1-dim LI velocity distribution for a relativis¬ 
tic gas. In Fig. [J] we plot the distribution ([ 8 ]) for a fixed mass and 
several decreasing parameters (3. 

The LI velocity distribution changes the concavity at zero veloc¬ 
ity for the particular value /3 = 1 which separates highly relativis¬ 
tic domains from domains that contain the Maxwell-Boltzmann limit. 

Above /? = 1 the zero velocity concavity is negative; below (3 = 1 
the concavity is positive characterizing an extreme relativistic limit 

Hn fact any inertial frame can be considered with respect to which the velocities are 
constructed. The results of the numerical simulations and the related histograms will not 
change as can be verified. 
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Figure 4: Plot of the ID LI velocity distribution ([8]) for a fixed mass and 
decreasing values of the adimensional parameter 13 (successively dashed, 
dashdotted and continuous curves, cf. text). The distribution is zero for 
> 1, as expected. 


for /3 <C 1. On the opposite side for /3 :§> 1 we have approach the 
MB regime. For the Jiitnner distribution Q with d = 1 the concavity 
at zero velocity changes for mo/3j = 3 showing that a high relativis¬ 
tic regime is attained at a lower temperature, as compared to the LI 
distribution (l 8 l) . 

We should finally remark that the variable a is in fact a particular 
case of the LI distance measure in the 3-dim Lobachevsky velocity 
space, which is the space of velocities in the special theory of relativity, 
as discussed in [5]. 

3 Discussion on how to construct his¬ 
tograms of velocities for a relativistic 
gas 

In recent past years a large number of numerical simulations of a rel¬ 
ativistic gas has appeared in the literature piisKiiKiaiiaiiT]. In 
these simulations the authors have studied some proposals of distribu¬ 
tion of velocities for a relativistic gas, like the Jiittner distribution [TO] 
or the Jiittner modified distribution [ 6 ]. The overall conclusion of all 
of these papers is that the Jiittner distribution matches accurately the 
data obtained by relativistic molecular dynamic simulations, convinc¬ 
ing the majority of the scientists that Jiittner is the correct relativistic 
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distribution. In the two following subsections we discuss some issues 
connected to the Lorentz invariance in these simulations and we dis¬ 
tinguish approaches to construct the histograms of velocities to be 
compared with the theoretical distributions. 

3.1 Histograms with constant velocity bins 

In the molecular dynamical simulations used previously in the above 
cited papers, the histograms of velocities have been constructed using 
the event driven simulation method [3 [8]. After starting the simula¬ 
tions with the particles of the gas having random position and veloc¬ 
ities and waiting the equilibration time, they collect the velocities of 
the particles at a precise moment and store them in a pool. The pro¬ 
cedure is repeated many times and when the pool of velocities has a 
large number of data a histogram is constructed having in the abscissa 
equal size bins of difference of velocities. In our simulations here we 
have 100 samples, each with N = 5, 000 particles and equilibration 
time equal to 100 x N. 

Let us discuss the construction of these histograms in one dimen¬ 
sion. The extension for two or three dimensions is straightforward. In 
one dimension we divide - in the abscissa - the domain of the values 
of the velocities obtained, say [—c, c], in n intervals (bins) of equal size 
(n = 60). In the ordinate of the j-th bin (j = 1, • • • ,n) we put the 
number of all particles having, at the moment of the measurements, 
velocities comprised between the limits of the bin. The histogram thus 
obtained, after normalized, is compared with the theoretical relativis¬ 
tic distributions proposed in the literature, as the Jiittner distribution, 
that in d-dimension is given as 

fj{v, m, 13j) = 7(v)2+'^ exp[-/3jmo7(v)] , (9) 

where mo is the rest mass of the particles and f3j = c^/kbT, where 
Kb is the Boltzmann constant. The one dimensional Jiittner distribu¬ 
tion is obtained by simply taking d = 1 in Eq. (j9|!FI. Adopting this 
procedure, all the numerical simulations in the literature referred to 
here match accurately with the Jiittner distribution (cf. for instance 
[ 6 ] in 1 -dim, [m [l 2 ] in 2 -dim and m in 3-dim). As a consequence 

^In the l-dim case it is necessary that a fraction of the N particles considered in the 
simulations (let us say N/2) have distinct masses in order that the distribution of initial 
velocities do not remain unaltered and that the system undergo equilibrationj^ [6]. 


10 



of these results the present well-established position adopted by most 
of the scientists of the statistical mechanics community is that the 
correct relativistic distribution is the Jiittner distribution. A striking 
illustration of this can be seen in (Figure 1), [11] (Figure 1) and 
[13] (Figure 1). 

In our simulations we have also performed a 1-dim relativistic dy¬ 
namical molecular simulation and the comparison of the histogram 
(constructed with constant velocity bins) with the Jiittner distribu¬ 
tion is displayed in Figures [D and [2|(/e/t) for heavy and light particles. 
As it can be seen, the agreement is also quite accurate. 

However we consider that two points are not satisfactory within 
this approach. The first one is that the theoretical Jiittner distribu¬ 
tion is not a LI distribution. In fact not only the relativistic energy, 
present in the argument of the exponential of Eq. ([9|), is not a LI 
quantity but also the correct 7 (v) factor needed to guarantee Lorentz 
invariance in 3D is 7 ^(v) (connected to the invariant element of vol¬ 
ume of the Lobachevsky space of relativstic velocities) and not 7 ^(v) 
as it appears in the Jiittner distribution. The second point is that, 
in the construction of the histograms, in the abscissa, where the //-th 
bin (which initially is constant for any /i and equal to A = 2c/n, for 
n bins) corresponds to particles having velocities going from to 
v^, is not LI, since the difference = v^ — (the size of the /z-th 
bin) changes under a Lorentz transformation, i.e., all the bins in the 
abscissa change differently under a Lorentz transformation. There¬ 
fore, the diagram is not LI as well, presenting a different form in each 
reference frame. Consequently, we have the following (uncomfortable) 
situation: a histogram that is not LI matches well a distribution of 
velocities (Jiittner) that is not as well LI. This certainly is not a satis¬ 
factory theoretical framework. What we have to search is to construct 
a LI histogram of velocities and then try to find a LI distribution of 
velocities that matches it. This is the aim of our work. 

We would like to remark that in a non-relativistic gas, whose par¬ 
ticles obey the Maxwell-Boltzmann distribution, the whole system 
should be invariant with respect to a Galilean transformation of ve¬ 
locities. In a one dimension gas, for simplicity, let us consider the 
difference Vi of velocities of the z-th particle (i = I, • • • , A) with the 
velocity of the box in a Galilean inertial frame, Vi = Vi — Ubox- In 
the abscissa let us divide the domain of the values of the variable Vi 
obtained in the simulations (say [—r'miru i^max] C (— 00 , 00 )) in n inter¬ 
vals (bins) of equal size where the /z-th bin (/z = I, 2, • • • n) is written 
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‘^{1 = (= A for all /r), with uq = -Umin and = /iA + z/q- 

The size of each bin is equal to (z/min + i^max)ln, where i^min is the 
maximum velocity obtained in direction —x and Vmax is the maxi¬ 
mum velocity obtained in direction +x. We remark that A^, for any 
/i, remains constant after a Galilean transformation. If in the ordi¬ 
nate we put in the /r-th bin the number of particles whose associated 
value of Ui = Vi — Ubox is comprised between z/^_i and the ordinate 
also does not change after a Galilean transformation. This implies 
that the whole histogram thus constructed is Galilean invariant and 
is well-fitted with a Gaussian centered at zero in any Galilean frame. 
This is the framework we want to reproduce in the relativistic case. 

We also note that if, in the ordinate, we put in the /r-th bin the 
number of particles having velocities comprised between and 
(not the difference of velocity of the particle and the velocity of the 
box as previously), after a Galilean transformation the histogram fits 
well with a Gaussian centered at a nonzero velocity, that is the veloc¬ 
ity associated with the Galilean transformation. Therefore, the best 
way to construct a histogram that is Galilean invariant is considering 
the difference of velocities between the particles and the box. The 
histogram so constructed fits with a Gaussian centered at zero in any 
Galilean frame. 

3.2 LI histograms: constant rapidity bins 

In the relativistic molecular dynamics, we also use the event driven 
simulation method with the relativistic scattering rules among the 
particles. In order to construct a histogram that is relativistic, or 
Lorentz, invariant, following the scheme of the non-relativistic case 
(that is Galilean invariant), we want that nor the abscissa neither 
the ordinate of the histogram change after a Lorentz transformation. 
Glearly, the bins constituted by the difference of velocities A^ = — 

are not LI. It is therefore crucial to use in the abscissa a quantity 
that is LI. To do this, let us first consider the relativistic difference d* 
between the velocity of the z-particle (vj) and the velocity of the box 
(tJfeox), given by (cf. dS])) 

2 _ (VZ - Vbox)^ 

(1 - Vi • Vfeoj,/c2)2 ■ 
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We now consider the rapidity, Sj, of this difference, concerning the 
particle i, that is a LI quantity, 

^ = tanh~^(p7c) ■ (11) 

We then propose to construct the abscissa of the histogram using the 
values of the rapidity, that now goes from minus to plus infinity. In 
the simulations we take the minimum and maximum numerical values 
obtained for the rapidity in all samples and divide this interval in 
the abscissa in equal rapidity-bins, where in the /x-th rapidity-bin we 
put all particles having values of {sj} going from the inferior and 
superior limits of the bin. The normalized histogram constructed in 
this way is invariant by a Lorentz transformation, as expected, but 
it is a histogram of the rapidity. In order to get a histogram of the 
velocities, or of the difference of velocities between the particle and the 
box, we use the inverse transformations of Eqs. (jlip and (IlOp (once we 
know the velocity of the box) on the abscissa. These transformations 
will change the relative size of the bins, since now the bins are no 
longer equally spaced in a velocity scale, as they are in a rapidity 
scale. However, even if we now represent the rapidity histogram in a 
velocity-scale, this histogram is still LI. 

In Figure [3] we show the histogram in a scale of difference of veloci¬ 
ties for N = 2, 500 light particles (the simulations have also N = 2, 500 
heavy particles) obtained from the histogram of the rapidities for the 
light particles. The agreement of the numerical data with the best fit 
of the theoretical LI velocity distribution ([8]) is excellent (continuous 
line), with a normalized rms between the curve and the points ~ 0.008. 
We also represent in Figure [3] the best-fit of the Jiittner distribution 
(dashed line) with the data of the histogram, where the parameter 
mo/3j has been adjusted. Clearly the Jiittner distribution does not 
fit the data, in contrast to the case where the histograms are con¬ 
structed, from the beginning, with constant bins of velocities. This 
indicates that the Jiittner distribution cannot be considered as the 
relativistic partner of the Maxwell-Boltzmann distribution. We would 
like to remark that in our simulations we obtained, with the method 
of constant bins of velocities, a good agreement with the results of [6], 
matching the Juttner distribution. Our results were obtained from 
100 simulations with N = 5 ,000 particles, with a simulation time of 
T = 100 X N and with a fixed number 60 of event points for the equi¬ 
librium conhguration. The best fit of the Juttner distributions yields 
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values of the parameter mQj3j with a relative difference of only 3% 
for heavy and light particles (cf. Figs. [T](2] (left)). However, larger 
differences of the values of (3 appear in the LI distributions and we are 
performing further simulations to completely clarify this point. 


4 Conclusions 

We have examined here the full relativistic molecular dynamics of a 
1-dim gas, by basically reproducing the numerical simulations made in 
a landmark paper on the subject [6], where the authors established nu¬ 
merically that the Jiittner velocity distribution function is the correct 
generalization of Maxwell’s velocity distribution in special relativity. 
We were led to look into this problem again since we considered such 
conclusion not satisfactory mainly due to the fact that the Jiittner dis¬ 
tribution is not LI. We therefore approached this problem using the 
same numerical simulations but adopting two distinct procedures on 
constructing the histograms with the data obtained. The first one cor¬ 
responds to the analysis of [6], where the histograms are constructed 
with equal size bins in the velocity variables and is - by construction 
- not LI. Using this procedure we were able to reproduce the results 
of [6], where the equilibrium distribution of the velocities resulting 
from the relativistic simulations fits accurately the Jiittner distribu¬ 
tion. The second procedure, which is a new contribution of this paper, 
adopts histograms based on the rapidity, a variable which is LI and is 
associated with the relative velocity of the particles: the histograms 
are constructed with equal size bins in the rapidity and are LI in the 
same way that the histograms of Maxwell-Boltzmann distributions, 
constructed with equal size bins in the velocity, are Galilean invariant. 
From these rapidity histograms we get a histogram of the relative ve¬ 
locities, or of the difference of velocities between the particle and the 
box, by using the inverse transformations of Eqs. (Ill|) and (jlOp (once 
we know the velocity of the box) on the abscissa. These transforma¬ 
tions will change the relative size of the bins, since now the bins are 
no longer equally spaced in a velocity scale, as they are in a rapid¬ 
ity scale. However, even if we now represent the rapidity histogram 
in a velocity-scale, this histogram is still LI. We then show that the 
equilibrium distribution of velocities is accurately fitted by the LI ve¬ 
locity distribution ([8|) proposed in this paper but not by the Jiittner 
distribution. In this way we believe that the symmetric connection 
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between the Maxwell-Boltzmann velocity distribution and the LI ve¬ 
locity distribution Q is established, the first Galilean invariant and 
the second LI, by using the variable relative velocity with respect to a 
given inertial frame, provided that the evaluation of relative velocities 
involves respectively the “arithmetic addition” or the “Lobatchevsky 
addition”. 
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